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Abstract Higher variability in rainfall and river discharge 
could be of major importance in landslide generation in the 
north-western Argentine Andes. Annual layered (varved) de- 
posits of a landslide dammed lake in the Santa Maria Basin 
(26°S, 66°W) with an age of 30,000 14 C years provide an 
archive of precipitation variability during this time. The com- 
parison of these data with present-day rainfall observations 
tests the hypothesis that increased rainfall variability played 
a major role in landslide generation. A potential cause of such 
variability is the El Nino/ Southern Oscillation (ENSO). The 
causal link between ENSO and local rainfall is quantified by 
using a new method of nonlinear data analysis, the quanti- 
tative analysis of cross recurrence plots (CRP). This method 
seeks similarities in the dynamics of two different processes, 
such as an ocean-atmosphere oscillation and local rainfall. 
Our analysis reveals significant similarities in the statistics of 
both modern and palaeo-precipitation data. The similarities 
in the data suggest that an ENSO-like influence on local rain- 
fall was present at around 30,000 14 C years ago. Increased 
rainfall, which was inferred from a lake balance modeling in 
a previous study, together with ENSO-like cyclicities could 
help to explain the clustering of landslides at around 30,000 
14 C years ago. 



1 Introduction 

Climate is a major influential factor for m ass movements in 
high m ountain regions. Increased humid ity JPethier and Reneaul 
19961) or increased variability in rainfall dGrosiean et al.ll997t 
Keeferetal 119981) can reduce thresholds for catastrophic land- 
sliding. In order to estimate the influence of climate in a given 
region, the climatic conditions during episodes with enhanced 
landsliding are compared with periods without important mass 
movements. The precise definition of climate scenarios at 
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times of high rock avalanche activity helps us to define thresh- 
old values for increased landsliding. 

About 30,000 14 C years ago, multiple large rock 
avalanches with volumes in excess of 10 6 m 3 occurred in 
the arid to semiarid intra-andean basins of north-western Ar- 
gentina dStrecker and MarretJl999UHermanns and StreckeJl999t 
iTrauth and Strecketil999l) . A potential mechanism that could 
have caused this enhanced landsliding in such an environment 
is climate change. Increased humidity and/or higher inter- 
and intraannual rainfall variability results in higher river dis- 
charge and erosion in narrow valleys and therefore increased 
destabilization of mountain fronts. 

The climatic conditions in NW Argentina are not well 
known for the period at around 30,000 14 C years ago. Marine 
and terrestrial records from tropical and subtropical South 

America indeed suggest more humid conditions (the Minchin 

period between 40,000 and 25,000 14 C years ago, e. ^. |van der Harnmen and A 

ll994tlLedru et al.ll996tlGodfrev et al.ll997llTurca et alll997l) 

and a strong El Nino/ South ern Oscillation (ENSO) (e. g. lOberhansli et"afl 

Il990t iBeaufort et aljEool . Various modeling studies have 

shown an impact of orb ital forcing on ENSO and its weaken - 

ing during the ice ages dClement et all 1 9991 Eiuet al.l2000l) . 

Thus El Nino events may be rare around 30, 000 14 C years ago 

dClement et al.ll999l) . lBeaufort et alJfcOOll) . however, inferred 

from Coccolithophores produ ction a significant occ urence of 

the ENSO for this period, and lTudhope et al.ld2.00ll> based on 

the analysis of annually banded corals concluded that ENSO 

has been a persistent component of the climate system over 

the past 130 ka. 

However, the local signal of the climatic changes in NW 
Argentina is still not well defined. Laminated sediments from 
a former landslide dammed lake in the Santa Maria Basin 
(NW Argentina, 26°S 66°W) contain valuable information 
about the environme ntal conditions for the period around 
30,00 I4 C years ago dTrauth and StreckeJl999tlTrauth et alJ 
l2000h . Hydrologic modeling of this palaeo-lake indeed in- 
dicates significantly wetter conditions during this time com- 
pared to the present (ar ound 10 to 15 % higher precipitation, 
iBookhagen et alj|200ll) . Linear spectral analysis of palaeo- 
precipitation data derived from annual layered (varved) lake 
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sediments also suggest an ENSO-like influence on rainfall leads to a typically dry winter climate (less than 50 mm per 
and consequently increase d interannual rainfall variability in m onth). 

river d ischarge and erosion dTrauth and Streckerll999HTrauth et all On interannual time scales, summer precipitation in the 

Central Andes, is primarily related to changes in meridional 
baroclinicity between tropical and subtropical latitudes, 
which in turn is a response to sea surface temperature anoma 



2000). However, the results of such linear methods are often 
ambiguous and not appropriate, since natural processes are 
complex, exhibit nonstationarities and are mostly recorded as 
short and noisy time series. In fact, data gained from sed- 
imentation processes (as colour data) are nonstationary by 
their origin and the relationship between climatic forcing and 
rainfall is not expected to be linear. Linear methods are usu- 
ally unsuitable to investigate natural complex data. In addi- 
tion, these approaches do not provide any information about 
a change in the climate dynamics, e. g. the sign of precipita- 
tion changes related to ENSO-like oscillations. 

The aim of our work is to test the hypothesis that an en- 
hanced ENSO-like influence on local rainfall caused a tempo- 
ral landslide cluster 30,000 14 C years ago. For this purpose, 
we first try to identify ENSO-like patterns in present-day pre- 
cipitation data and infer causal links between this ocean-atmos 
phere oscillation and local rainfall. Secondly, we search for 
similar influences in palaeo-precipitation data reconstructed 
from 30,000 14 C year old lake sediments. This comparison 
is carried out using a new method of nonlinear data analysis, 
the cross recurrence plots (CRP), whic h can be applied to 
short and nonstationary complex data jMarwan and Kurthsl 
120021) . This procedure traces similarities and differences in 
several measures of complexity in both modern and past rain- 
fall data. Significant occurrences of the ENSO-rainfall tele- 
connection together with increased rainfall could help to ex- 
plain enhanced landsliding 30,000 14 C years ago in NW Ar- 
gentina where no major mass movements occur today. 



2 Present-day Climatic Conditions 

Summertime climate and atmospheric circulation over NW 
Argentina is large ly governed by the S outh American mon- 
soon system (e. g. Z hou anc lLaulll998h . featuring heavy pre- 
cipitation, an upper-air anticyclone (Bolivian High) and a low- 
level trough (Chaco low). Approximately 80% of the an- 
nual precipitation am ount falls within the sum mer months 
November -February dBianchi and Yane3ll992l) . associated 
with southward moisture transport t o the east of the Andes 
throu gh the Andean low-level jet (e. g. lNogues-Paegele and Mol 
Il997n . The intra-andean basins and valleys, separated from 
this low-level moisture flux through the north-south running 
eastern Andean ridge, are arid and receive less than 
200 mmyr 1 , whereas the regions to the east of this pro- 2M 
graph ic barrier receive more than 1 500 mm yr~ 1 dBianchi and YanezT 



lies in the tropical Pacific (e. g. l Vuille et al]200otlGarreaud and Aceituno| 
l200lt iGarreaud et al J lin pr ess). The study region therefore 
shows a significant relationship with ENSO, featuring a weak- 
ened westerly flow with a significantly enhanced easterly mois- 
ture transport during La Nina summers and strengthened west- 
erly flow with a significantly subdued easterly moisture trans- 
port during El Nino summers. As a result, the rainy season 
is much more active during La Nina episodes and less active 
during El Nino episodes. These 

ENSO related atmospheric circulation anomalies are also ev- 
ident in radiosonde data to the east of the Central Andes over 
NW Argentina (Salta), featuring enhanced southeasterly 
(northwesterly) flow and increased (decreased) specific hu- 
midity levels in the lowe r and mid-tr oposphere during La 
Nina (El Nino) summers JVuillal 1 9 99 > . The notion that this 
ENSO influence indeed extends beyond just the Central An- 
des is further supported by several recent studies, which in- 
dicate that precipitation anomalies in this part of the Andes 
tend to coincide with anomalies of the same sign over SE 
Bolivia and NW Argentina (e. g. [A ceituno and MontecinosI 
ll997tlGarreaudll999l) . lBianchi and Yafiezl Jl992l) also report 
a weak but significant tendency toward less rain during El 
Nino ye ars, based on a high -density network of 380 weather 
stations. iTrauth e t al. ( 2000) provide a detailed statistical anal- 
ysis on the same data set showing that this tendency is very 
obvious but spatially and temporally highly variable. As indi- 
cated in Figure^for the representative El Nino event 1965/66, 
precipitation can be decreased up to 80% with respect to 
the long-term average in rainfall, with a more significant re- 
duction in the northern part of the study area. A composite 
analysis of the monthly summer precipitation (DJFM) differ- 
ence between El Nino and La Nina summers between 1979 
and 1999 based on CM AP satellite-derived precipitation data 
dXie andArkinll997l) confirms this notion (Hig. |2j. Although 
weak, the tendency toward less precipitation during El Nino 
and more precipitation during La Nina episodes is evident 
even in this low-resolution gridded data. Similar to the pat- 
tern in Figure the ENSO signal is reversed a few degrees 
further south, where summer precipitation becomes less dom- 
inant. 



3 Methods 



I1992I) . 

Due to the seasonal change in the tropospheric tempera- 
ture gradient between low and mid-latitudes, the subtropical 
westerly jet extends further north during the winter months, 
reaching its northernmost position around 27°S. The resulting 
wintertime mean westerly flow, which prevails over the study 
region in the mid- and upper troposphere, hinders regional 
moisture transport over the eastern slopes of the Andes and 



It is difficult to compare rainfall proxies from 30,000 14 C 
year old lake sediments with present precipitation data. The 
process recording weather and climate in p alaeo-archives is 
complex and s o far not very well understood dSaltzmanll990t 
iBradlevll 19991) . Because of various signal distortions in both 
time and frequency domain, the use of linear data analysis 
methods reaches its limits. The complexity of natural pro- 
cesses suggests the application of nonlinear methods instead 
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Fig. 1 Study area in the Santa Maria Basin with the locality of an- 
nual layered lake deposits in the locality El Paso, the relative precip- 
itation anomaly during the El Nino 1965/66 compared to mean an- 
nual preci pitation (annual precipitat ion as a mean from July to June; 
data from lBianchi and Yafiej[l993) and the prevailing w ind direc- 
tions during January (black arrows; wind directions from lProhaskal 
Il97d) . 
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Fig. 2 Difference in precipitation (in mm day -1 ) between El 
Nino and La Nina periods (El Nino - La Nina) for December 
— March based on C PC merged analysis of precipitation (CMAP, 
Xie and Arkii ] | 19971) between 1979 and 1999. Contour interval is 
—4, —2, — 1, —0.5, 0.5, 1, 2, 4 mm day -1 , regions with difference 
> 0.5 (< —0.5) shaded in light (dark) gray. An El Nino (La Nina) 
event is defined as a phase of at least six consecutive months in 
which the 5-month running mean of SSTA in the NIN03.4 region 
exceeds (or remains below) 0.5°C (— 0.5°C). The reference period 
for the SSTA is 1961-90. 



for the analysis and comparison of such complex processes. 
Most of the nonlinear standard techniques, such as fractal di- 
mensions or Lyapunov exponents, ca nnot be estimated for 
such data dKantz and Schreiberlll997|) . Therefore, we have 
tried to quantify cross recurrence plots of present-day and 
palaeo-data. This reveals a suite of complexity measurements 
which give hints to identify similar patterns in present-day 
and palaeo-data. This comparison first tests the hypothesis 
that the signal extracted from the lake sediments is an ap- 
propriate measure for palaeo-precipitation. Secondly, it helps 
to reconstruct the variability in annual rainfall as compared 
to the present. Both results can then be used to test the hy- 
pothesis that increased interannual variability in climate can 
explain enhanced landslide activity 30,000 14 C years ago. 



3.1 Cross recurrence plots ( CRP) 

An important aspect of climate changes involves nonlinear 
interactions among many compo nents of the earth's environ- 
mental system (e. g. |PalmeJl999h . These components include 
the oceans, land, lakes and continental ice sheets, and in- 
volve physical, biological, and chemical processes. Many of 
the techniques used to diagnose climatic variability such as 
Fourier analysis, empirical orthogonal functions or singular 
value decomposition are formulated using methods taken from 
linear analysis. However, while using these techniques it is 
difficult to analyze the nonlinear character of the earth's cli- 
mate system. In the last two decades, a great variety of nonlin- 
ear techniques have been developed to analyse data of com- 
plex processes. Most popular are methods to estimate fractal 
dimensions or Lyapunov exponents (e. g.lMandelbrotll 19821 
IWolf et al.ll985tlKurths and HerzeJl987tlKantz and Schreibei 
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Il997h . However, a number of pitfalls are possible due to the 
uncritical use of these methods on natural data, which are 
typically nonstationary and noisy. Furthermore, we cannot 
expect low dimensions in highly complex natural processes. 
Therefore, we have modified and applied the nonlinear data 
analysis method of cross recurrence plots, which was recently 
introduced by the extensio n of recurrence plots JZbilut et alJ 
ll998llMarwan et all2002l) . for detecting similarities and dif- 
ferences in the ENSO influence in present-day and past rain- 
fall data. In order to quantify such similarities by using CRPs, 
new measures of complexity were introduced. Measures of 
complexity were developed in order to quantify the complex- 
ity of processes; the simplest measure is the entropy, which 
distinguishs between noisy and periodically processes. Here 
we use an approach which uses the geometrical structures 
which are contained in CRPs. This new technique is partic- 
ularly efficient for the analysis of rather nonstationary, short 
and noisy data and was successfully applie d to prototypical 



The CRP is a two-dimensional N x N array of points 
where N is the number of embedding vectors obtained from 
the delay coordinates of the input signal. The values of the 
CRP are one (black points) if trajectories lie close to each 
other (recurrence points), whereas values of zero (white 
points) document rather large distances between two trajec- 
tories. From the occurrence of lines in the CRP parallel to 
the diagonal in the recurrence plot it can be seen how fast 
neighbouring trajectories diverge in the phase space. Recur- 
rent data in a system would create diagonal lines in a dis- 
tance t from the main diagonal in such a plot comparing both 
phase-space embeddings with respect to the time delay t . It is 
important to note that an additional CRP with opposite signed 
second time series CR^- = ®(s — ||x,- + y / ||) allows to distin- 
guish positive and negative relations. 

Visual inspection of CRPs already reveals valuable in- 
formation about the relationship between two complex pro- 
cesses. However, a better understanding of causal links be- 



mode l systems with nonlinear interrelations JMarwan and Kurtha tween both processes demands a more quantitative exami- 
nation of the CRPs. Therefore, we intr oduce the followin 
two s tatistical measures of complexity ( iMarwan and Kurthi 
120021) : 

the recurrence rate, 



2002). Thus, the CRP is an appropriate method for time series 
analysis of climate and palaeo-climate data. 

The basic idea of this approach is to compare the dynam- 
ics of two processes which are both recorded i n a single tim e 
series. Following Taken's embedding theorem llTakensll98ll) . 
the dynamics of a process with m state parameters (i.e. m 
differential equations), which is, however, measured by only 
one time series u(t) — Uj with length N and a sampling time 
At (i.e. t — iAt), can be appropriately presented in its recon- 
structed phase space of a dimension m (theoretically when 
m > 2m + 1). Such a reconstruction can be formed by using 
the time delay method (embedding), where for each compo- 
nent of the state vector a value of the time series after a pre- 
defined delay T (time delay) is choosen: 
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The dimension m of such a reconstructed state or phase 
space is called embedding dimension. The time evolution of 
the state vectors forms a trajectory x,, which runs through all 
possible states at time t = iAt and, thus, present the dynamics 
of the process. 

The similarity in the behaviour of both processes in this 
reconstructed phase space can be examined by using the CRP, 
which visualizes the distance between segments of their phase 



where Z mi „ is a predefined minimal length of a diagonal line 
segment, P (I, t ) is the histogram of the diagonal line lengths 
in CR ± at a distance t from the main diagonal (i. e. the time 
delay t between the two phase space vectors) and t = iAt. The 
RR(t) is the density of adjacent states, i. e. of the recurrence 
points in a CRP diagonal. RR(t ) therefore measures the prob- 
ability of similar states in both processes after a delay time t . 
High densities of recurrence points in a diagonal cause high 
values of RR, which is typical for processes with a similar 
behaviour in the phase space. 

Strongly fluctuating data cause short or absent diagonals 
in the CRP, whereas data from deterministic processes pro- 
duce long er diagonals. If two deterministic processes have the 
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space trajectories x, and y, of the embedded time series JMarwan aadiifairthaimilar phase-space behaviour, i. e. the phase-space 

trajectory reaches the same regions of the phase space dur- 
ing certain times, then the number of longer diagonals will 
increase and the amount of short diagonals decrease. The av- 
erage diagonal length L measures the epoch length (i. e. the 
time span) of significant similarities in the behaviour of both 
processes. The higher the coincidence of both processes, the 
larger the length of these diagonals. 

Consequently, high values of RR and L correspond to fre- 
quent and longer periods of similar behaviour of the pro- 
cesses as recorded in the time series data. Therefore, these 
parameters are appropriate quantitative measures for the sim- 
ilarities between both processes. However, extrema at longer 



where e is a predefined cut-off distance, || • || is the norm 
(e.g. the Euclidean norm) and ®(x) is the Heaviside func- 
tion. Depending on the type of the application, e can have a 
fixed value or can vary for each ; in such a way that a pre- 
de fined number of neighbours occur within a certain radius 
e dEckmann et al.lll987HMarwan et~al]l2002l) . This results in 
a constant density of recurrence points in each column of the 
CRP and is particularly useful in the analysis of complex pro- 
cesses with differences in the variability of the amplitudes. 
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delays t do not necessarily correspond with high correlations. 
Future work will concentrate on the theoretical and more de- 
tailed investigation of the interrelations between the struc- 
tures in CRPs. 

We have proposed a statistical evaluation of the quantita- 
tive measures of the CRP with an ensemble of a large amount 
of surrogate data. 

The assumption for the surrogate data is that the consid- 
ered processes are linearly independent and do not have any 
similar dynamics. These surrogates should reveal some fea- 
tures like in our original data but also features caused by the 
randomness of a possible correlation (stochastic processes). 
Linear cor related noise is a paradigm atic example for such 
processes jKantz and Schreiberlll997l) . We calculate a surro- 
gate time series based on this class of processes with the fol- 
lowing recursive function, a autoregressive process of order 
Pi 

p 

x n — cipX n _p -\- b^ n: 

i=\ 

where t, is white noise and a, are coefficients which de- 
termine the auto-correlation of the system and allow to adapt 
this stochastic system to our natural processes. We fit the 
model to the precipitation series of the station Tucuman. Then 
we perform the CRP analysis using the SOI data and the 
ensemble of, e.g. 10,000 realizations of precipitation series 
produced by the AR model. Using the distributions of the RR 
and L measures we can estimate their empirical confidence 
bounds (we will use the 2c-bounds which approximately cor- 
respond with the 95% confidence level). 

With these confidence bounds we can evaluate the ob- 
tained measures of CRP and the relations of the natural pro- 
cesses. Since the surrogates are from a stationary system and 
the natural data are nonstationary, we have further applied 
this kind of evaluation to more stationary segments in the nat- 
ural data. We got the same results. This kind of surrogates is a 
special realization, which is prototypical for linear stochastic 
processes, and there are a lot of other possibilities to construct 
surrogates. 



3.2 Comparison of modern and palaeo-precipitation 
variability 

In order to test the new method on precipitation data, we first 
compute the CRP for rainfall stations with well established 
and clear ENSO influence. We use monthly precipitation data 
from the cities of Buenos Aires (BAI) and Caracas (CAR) 
from the WMO data set JHoffmannlll975l) . For the assess- 
ment of the modern ENSO influence on local rainfall in NW 
Argentina, we analyze monthly precipitation data from three 
stations: San Salvador de Jujuy (JUY), Salta (SAL) and San 
Miguel de Tucuman (TUC; Figs. ^ an d|3]B). These stations 
in the capitals of the provinces Jujuy, Salta and Tucuman 
provide the longest time series from this area and are lo- 
cated on a north-south transect. Moreover, these locations are 
influenced by different local winds; Jujuy and Salta mainly 
receive north-easterly and easterly moisture-bearing winds 



during the summer rainy season, whereas Tucuma n is char- 
acteri zed by southerly and south-westerly winds dProhaskal 
11974 . The Southern Oscillation Index (SOI) is used as a 
measure for the ENSO variability between the years 1884 
and 1990 (Fig. |3]A, based on COADS data). This index is 
the normalized difference between the sea level air pressure 
in Tahiti and Darwin. Extreme negative values represent El 
Nino e vents and extreme positive valu es represent La Nina 
events (IRopelewski and Halpernll987ft . In our analysis, we 
use monthly data, i. e. twelve data points per year in order 
to avoid loosing valuable intraannual information. Moreover, 
longer data vectors improve the significance of the CRP mea- 
sures, thus the use of annual data would significantly reduce 
the value of our results. The potential distortion of the final 
result by differences in the causal linkage between ENSO and 
rainfall over the year, i. e. between dry and wet season, is be- 
lieved to be of minor importance since the absolute values 
of precipitation during the dry season are low and hence the 
contribution to the analysis is small. The rainfall variability 
during the dry season is not significantly different from white 
noise and disappears after low -pass filtering preceding the ac- 
tual time series analysis. 
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Fig. 3 Smoothed and o-normalized time series of the South- 
ern Oscillation Index (A), monthly precipitation data of Salta 
(B) and its smoothed and o-normalized time series (C). SOI 
based on COADS data from the NOAA Live Access Server 
I http://ferret.wrc.noaa.gov l. 
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The CRP analysis of the present-day ENSO and precip- 
itation data reveals characteristic patterns that can now be 
traced in palaeo-precipitation data. The palaeo-precipitation 
variability is inferred from varved lake sediments sampled 
at the location El Paso (EP; 26.0°S, 65.8°W) in the Santa 
Maria Basin in NW Argentina (Figs. and |4j. These sedi- 
ments wer e deposited in a landslide dammed lake 30,000 14 C 
years ago JTrauth and Streckerl Il999t l iHermanns and Streckerl 
Il999h . Because of the internal structure of the deposits with 
intra-varved changing of diatom species and the cyclic recur- 
rence of pa ired diatomite and clastic l ayers, these laminations 
are varves JTrauth and Str ecker 1999). The annual cycle with 
wet summers and dry winters caused significant changes in 
the lake sedimentation. During the rainy season mainly ocher 
coloured silty sediments were deposited; during the subse- 
quent dry season a thin white layer consisting of the skele- 
tons of silica algae (diatoms) was deposited. Due to its white 
colour, the diatomaceous layers can be used to identify sin- 
gle years in these sediments. Recurring intense red coloura- 
tion of the silty part of the annual layers is sourced from re- 
worked older sediments which are eroded and deposited only 
during extreme rainfall events. Therefore, the intensity of red 
colour in the varved deposits can be in terpreted as a proxy 
for precipitation variatio n at El Paso site (Tr auth and Streckerl 
ll999HTrauth et al.l2000l) . The more intense the red colour the 
higher was the precipitation during the rainy season. 




Fig. 4 Detail of varved lake sediments from the El Paso site in the 
Santa Maria Basin with cyclic occurrence of dark red colourations 
recording more precipitation and sedi ment flux with ENSO-like pe- 
riodicities JTrauth and StreckeJl999T) . The overlayed curve shows a 
representative red colour intensity transect of the deposits. 



The colour intensity of a section of the sediments pro- 
file with a length of 160 varves was gained by scanning high 
quality photographs. After digital pre-processing, a time se- 
ries of red intensity values on a length scale was obtained. We 
transform these data to a time scale assuming an annual re- 
currence of the diatomaceous layers. Within single varves 12 
subannual data points are computed by logarithmic interpola- 
tion of the data taking into account the exponential decrease 
of the sedimentation rate during the annual cycle (Fig. [3J. 
The power spectrum estimate of the red colour intensity re- 
veals significant peaks within the ENSO freque ncy band of 2 
to 4 y ears, suggesting an ENSO-like influence JTrauth et alJ 
2000). Because of the nonstationarity of these data (the sedi- 
mentation process in a lake is not stationary, resulting in non- 
stationary proxy parameters for the in-lake processes; mean 
of the first half of the time series is 0.30, of the second half is 
—0.32; standard deviation of the first half of the time series 
is 1 . 13, of the second half is 0.71), linear correlation analysis 
is unsuitable. Therefore, we apply the CRP analysis to these 
data. 

250 1 1 1 1 1 1 
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Fig. 5 Red intensity values of the lake sediments of site EP160 on 
(A) a length scale and on (B) a time scale and after smoothing and 
normalization; the unit of raw data is one bit, the unit of transformed 
and smoothed data is the standard deviation a. 



4 Results of Nonlinear Data Analysis 

First, all time series are normalized and low-pass filtered us- 
ing a 7th-order Butterworth filter with a cutoff frequency of 
1/18 month -1 in order to remove the predominant annual cy- 
cle from the data (Figs. |5]C and [5}- Butterworth filters are 
from the infinite-duration impulse response type (IIR filters) 
and have a monotonically decreasing response with respect 
to frequency dElliotil987t> . 
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Fig. 6 Cross recurrence plot of SOI vs. precipitation data from the 
city of Salta (SAL). The x-axis shows the time along the phase space 
trajectory of the SOI and the y-axis that of SAL. Black points rep- 
resent the occurrence of similar states in both processes. Diagonal 
lines correspond with epochs of similar dynamics in both processes. 
The amount and length of these lines can be used as measures of the 
similarity of both processes. 



Next, the filtered rainfall data and the Southern Oscilla- 
tion Index (SOI) are embedded into a phase space using m = 3 
and x = 9. The method of nonlinear time series analysis us- 
ing delay time embedding relies on a choice of good delay 
time and the embedding dimension. Proper values for these 
parameters are determined using the m ethods of false near- 
est ne ighbours and mutual information JKantz and Schreiberl 
Il997l) . The quantitative analysis of cross recurrence plots is 
then applied to pairs of time series, local precipitation records 
and the Southern Oscillation Index (SOI). The CRPs are com- 
puted using a fixed amount of nearest neighbours with e = 
15%. Since the statistics of CRPs are sensitive to changes in 
the cutoff distance, we have run sensitivity tests in order to 
find the optimum value of e. The value of 15 % appears to be 
the best choice receiving robust and precise results. 

The CRPs of all pairs of time series show similar features. 
The significant similarities between CRPs obtained from 
modern (Fig. |6} and palaeo-precipitation data (Fig. Q indi- 
cate that the red colour intensity records from the varved lake 
sediment do reflect rainfall in NW Argentina. First we dis- 
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Fig. 7 Cross recurrence plot of SOI vs. the best matching section of 
palaeo-precipitation (EP160). Scaling as in Fig.© The x-axis shows 
the time along the phase space trajectory of the SOI and the y-axis 
that of EP160. 



cuss the CRP of Salta precipitation (data series SAL) vs. the 
Southern Oscillation Index (SOI) and the CRP of red colour 
intensity of varves (data series EP160) vs. SOI. The x-axis 
represents time along the phase space trajectory of the SOI, 
whereas the y-axis represents the time along the phase space 
trajectory of SAL or EP160, respectively. The CRP of SAL 
vs. SOI exhibits longer diagonal lines in two to four year in- 
tervals, which matches the same frequency band obtained by 
the power spectral analysis (Fig.^Ji. This indicates that some 
sequences of the phase space trajectory of the SOI recur in 
sequences of the phase space trajectory of SAL after relocat- 
ing by the time of two to four years. Vertical white bands in 
the CRP represent less frequent states in SOI, such as hori- 
zontal white bands suggest for SAL. The latter occurs with 
intervals of more than ten years. The CRP between EP160 
and SOI shows similar characteristics as the CRP described 
above (Fig.0. Longer diagonal lines have spacings of about 
two to four years. White bands occur at time scales of more 
than ten years. Some linkages in both CRPs are obvious by 
visual inspection. Next, the quantitative analysis of the CRPs 
is performed in order to study statistically these relations and 
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Fig. 8 RR and L measures of the cross recurrence plots between SOI 
and precipitation in Caracas (A, C) and Buenos Aires (B, D) with a 
well-established and clear ENSO influence. Extreme values reveal 
high similarity between the dynamics of the rainfall and the ENSO. 



to allocate the predefined causality patterns to certain locali- 
ties. 

In order to calculate the measures of complexity RR and 
L between the rainfall data and SOI, we used a delay time in 
the range between —12 and +22 months, i. e. these measures 
are determined in a small corridor above and below the main 
diagonal. We are interested in the extrema and in the time 
lag where they occur and we get the following results for the 
various pairs of records. From an ensemble of 10, 000 realiza- 
tions of a 5th-order AR-model we calculate the 2o-bounds of 
their distributions for RR and L. The coefficients for the AR- 
model are adapted to the Tucuman precipitation (we also used 
AR-models adapted to the rainfall data of the other stations, 
which revealed similar results). The order of the AR-model 
is determined with the Akaike's Information Criterion and a 
criteri on, which assesses whether the residues follow white 
noise JSchlittgen and Streitber gll999t) . 

The CRP measures between CAR and SOI reveal extreme 
positive values and between BAI and SOI extreme negative 
values, which reflect the strong influence of ENSO in these 
areas (Fig. [8). The parameter RR of the CRPs between TUC 
and SOI has small negative values, which do not exceed the 
2a-bounds, and does not show preferences for a distinct lag. 
The parameter L has also small values, but it has rather small 
maxima and minima at delays of —1, 4 and 8 months. These 
results indicate that the precipitation in Tucuman is not 
strongly influenced by ENSO. If there is a weak influence, 
the rainfall would increase during El Nino (Fig.|9]A, E). The 
analysis of JUY and SOI reveals clear positive values around 
a lag of zero and negative values after about 8-12 months, 



which suggests a significant link between Jujuy rainfall and 
ENSO (Fig.|9]C, G). The measures for the analysis SAL vs. 
SOI show smaller maxima for a delay of about zero and min- 
ima after a lag of about 8-12 months. We therefore infer a 
weaker linkage between Salta rainfall and ENSO (Fig.|9]B, F; 
the disrupted minima in the L parameter at around ten months 
is due to the short data length and a resulting nonstationarity 
in the CRP). The measures for SAL and JUY exceed the 2a- 
bounds. 

The 30,000 14 C year old precipitation data are not sim- 
ply comparable with present-day data, because there is no 
information available about how to synchronize the rainfall 
records with modern climate indices. Therefore we seek the 
time window in these data showing the highest coincidence in 
the dynamics using maximum values for RR and L as the key 
criterion. Although the observed coincidence is not very high, 
it yields the time section in the palaeo-precipitation record 
EP160 which can be best correlated with modern data. In 
our palaeo-data EP160 we find such a section represented 
by maximum and minima values for RR and L for delays of 
about zero and ten months, similar to those found for JUY 
and SAL (Fig.|9]D, H). The RR and L measures exceed also 
the 2o-bounds. 

To use the minima at lags around 8-12 months for cli- 
matological interpretations is difficult and might lead to er- 
roneous conclusions, but these characteristic patterns of pos- 
itive and negative interrelations can be used to compare the 
present-day and palaeo-data. The positive and negative inter- 
relations have the same time delay between 10 and 12 months 
in the present-day and the palaeo-data. 



5 Discussion 

We applied the method of cross recurrence plots (CRPs) to 
modern and palaeo-precipitation data in order to compare the 
magnitude and causes of rainfall variability in the NW Argen- 
tine Andes today and during the time of enhanced landslid- 
ing at around 30,000 14 C years ago. CRPs are able to look for 
nonlinear interrelations between two processes. The major re- 
sult from this analysis is the significant similarity between 
the complex dynamics of modern rainfall and the palaeo- 
precipitation as recorded in the red colour intensity record 
from the lake sediments in the location El Paso. The distances 
between longer diagonal lines in the CRP of both records are 
about two to four years, the approximate time of recurrence 
of extreme ENSO phases today. The first implication of this 
result is that the red colour intensity of the sediments is in- 
deed a good proxy for the rainfall intensity 30, 000 14 C ago. 
This r esult is in line with the observations of iTrauth et alJ 
fcoori suggesting an enhanced erosion of red-coloured clas- 
tic sediments during heavy rainfall events today whereas pre- 
cipitation usually only reaches the elevated areas with mainly 
greenish low-grade metamorphic rocks exposed. This effect 
causes predominant gre enish to buff-coloured clays deposited 
in the form er lake basin jTrauth and Streckerll999UTrauth et aD 
I2000l2003h . 
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Fig. 9 RR and L measures of the cross recurrence plots between SOI and precipitation in Tucuman (A, E), Salta (B, F), Jujuy (C, G) and 
palaeo-precipitation (D, H). Extreme values reveal high similarity between the dynamics of the rainfall and the ENSO. The dash-dotted lines 
are the empirical 2o-bounds from the distributions of an ensemble of data based on a 5th-order AR-model. 



Since our analysis of modern data reveals a strong rela- 
tion between local rainfall in the northern part of the study 
area (Jujuy) and ENSO, we interpret this similarity as an in- 
dication of a strong ENSO-like influence in the Santa Maria 
Basin at around 30,000 years. In contrast, there is no signif- 
icant linkage between the modern rainfall in Tucuman and 
ENSO. This result could indicate that ENSO does not influ- 
ence precipitation in the southern part of the study area or this 
influence is rather diffuse or changing in time. 

The CRP between the SOI and the rainfall in Jujuy and 
Salta reveals a positive relation without any large delay, i. e. 
the occurrence of an El Nino at the end of a year would cause 
decreased rainfall in the rainy season from November to Jan- 
uary and the occurrence of a La Niiia would cause increased 
rainfall during this time of the year. The opposite response 
after a delay of 8 - 12 months is not easy to interpret, be- 
cause we do not know which mechanism actually caused this 
linkage. The time span between the identified maxima and 
minima is about one year and could be explained by the fact 
that La Nina events often follow El Nino events. The smooth 
shape of the CRP measure curves are artefacts caused by 
low-pass filtering of the time series. The measures of CRP 
of Tucuman precipitation and SOI show non-significant val- 
ues without any characteristic delays. The analysis of varve 
data reveals a significant positive relation between SOI and 
palaeo-precipitation at the location El Paso. Similar to the 
modern situation, the CRP shows a significant negative rela- 
tion with SOI after a delay of about ten months. Both interre- 
lations are rather similar to those of ENSO-JUY and ENSO- 
SAL. 

The similarities between the time series of the modern 
rainfall data and the palaeo-precipitation record from the lake 
sediments suggests that an ENSO-like oscillation was active 



at around 30,000 14 C years ago (roughly corresponding to 
34,000 cal. years BP), which corresponds with t he results of 
the in vestigation of Coccolithophores production dBeaufort et al.l 
l200ll) . A younger landslide cluster in the same region at around 
5000 14 C (corresponding to 5800 cal. years BP) was also ex- 
plaine d by a stronger ENSO influe nce at that time JTraufh et al.l 

200C 1 palaeo-ENSO evi dence from lKeefer et al.ll998HSandweiss et ail 
200l UHau g etall200ll) . The spacing between both landslide 



clusters is around 28,000 years. Although two landslide clus- 
ters do not allow to infer a systematic recurrence of such 
events, we believe that there is some evidence that these events 
correspond to the periods of a strong ENSO-like variation as 



repor ted from deep-sea sediments off-shore Peru dOberhansli et al.l 
in the Indo-Pacific Ocean dBeaufort etalJl200lh and 



1199 

New Guinea corals jTudhope et al.ll200ll) . These long-term 
ENSO records suggest a mixed precession-glacial forcing on 
ENSO resulting in significant 23- and 30-kyr cyclicities, which 
confirms model results and recen tly inferred relations between 
ENSO variability and insola tion dClement et al.ll99"9llLiu etall 
l2000HRittenour et all200(il) . 



In the semiarid basins of the NW Argentine Andes, the 
ENSO-like variation could have caused significant fluctua- 
tions in local rainfall at around 30,000 14 C years similar to 
the modern conditions. Together with generally higher mois- 
ture levels as indicated by lake balance modeling results, this 
mechanism could help to explain enhanced landsliding at around 
30,000 and 5,000 14 C years ago in the semiarid basins of the 
Central Andes. 
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6 Conclusions 

The quantitative analysis of cross recurrence plots has re- 
vealed similarities in the evolution of the phase space tra- 
jectory of climate indices and present-day and past rainfall. 
In comparison to the usually less variable climate during ice 
ages, our result suggests an enhanced impact of ENSO-like 
conditions on local climate in the Santa Maria Basin 30,000 
14 C years ago associated with a strong inter- and intraannual 
variability of rainfall and an intensification of moisture trans- 
port. A more variable climate due to an enhanced ENSO-like 
impact could have raised the risk of landsliding in this re- 
gion and could help to explain enhanced landslide activity at 
around 30,000 14 C years ago. 
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